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Considering the defects of the previous work for estimating the anomalous production rates of — > 
T(lS)7r + 7r, r(2S)n + n- near the peak of the T(5S) resonance at yfs = 10.87 GeV [K.F. Chen et al. (Belle 
Collaboration), Phys. Rev. Lett. 100, 112001 (2008)], we suggest a new scenario where the contributions from 
the direct dipion transition and the final state interactions interfere to result in not only the anomalously large 
production rates, but also the lineshapes of the differential widths consistent with the experimental measurement 
when assuming the reactions are due to the dipion emission of Y(5S ). At the end, we raise a new puzzle that 
the predicted differential width dT(Y(5S) — > Y(2S)n + n~)/dcos9 has a discrepant trend from the data while 
other predictions are well in accord with the data. It should be further clarified by more accurate measurements 
carried by future experiments. 

PACS numbers: 14.40.Pq, 13.20.Gd 



Looking at the spectra of the bb system listed in the particle 
data book J2], there exist six bottomonia with J PC = 1 , 
which are, respectively, Y(nS) (n=l,2,3,4), Y( 10860), and 
T(l 1020). The first five resonances are orderly assigned as nS 
(n= 1, ■ ■ • , 5) bb states, whereas the extra one, T(l 1020), may 
be the 65 state. The estimate on their spectra in the potential 
model supports such assignments [2-4]. However, recently 
anomalous large rates of e + e~ — > T(lS)n + n~, Y{2S)n + n~ 
near the peak of the Y(5S) resonance at s/s — 10.87 GeV 
were observed by the Belle Collaboration to be larger than 
the dipion-transition rates between the lower members of the 
T family by 2 orders of magnitude. The Belle data are 
T(Y(5S) -» Y(lS)7r + 7r) = 0.59 + 0.04(stat) + 0.09(syst)MeV 
and r(Y(5S) -> Y(2S)7r + 7T) = 0.85 + 0.07(stat) ± 0.16(syst) 
MeV |5J. The Belle observation has stimulated theorists' ex- 
tensive interest in exploring the reason that results in such 
anomalous phenomena. 

There are two possibilities that may offer reasonable inter- 
pretations of the anomalous large rates. First, these anomalous 
production rates announced by Belle are from an exotic res- 
onance structure different from Y(10860). The second is that 
there may exist extra contributions that differ from the direct 
dipion emission Y( 10860) — > Y(15, 2S)n + n~. It is expected 
that with a careful analysis we may eventually identify the 
reasonable, or at least the dominant source of the large rates. 
Thus one not only obtains the branching ratios, but also needs 
to fit the lineshapes of the differential widths over the invariant 
mass of dipion and over the angular distribution cos 8. 

Along the first route, Ali et al. suggested a tetraquark 
interpretation of y 6 (10890) = [bq][bq] &[§ and ana- 
lyzed the Belle data |5] for the anomalous Y(lS)7r + 7r~ and 
Y(2S)n + 7T~ productions near the Y(5S) resonance. By fitting 
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the ji + n invariant mass spectrum and the cos 6 distributions 
for y 6 (10890) -> T(nS)T: + 7T- (n = 1,2) shown in Figs. 1 and2 
of Ref. 0], they claimed that the tetraquark interpretation can 
well describe the anomalous rates of the two-pion-production. 
In their scenario, there are non-resonant and resonant contri- 
butions interfering to result in the branching ratio and differ- 
ential widths. In that work 0Htl, a simple Lorentz structure 
Lagrangin was introduced to stand as the effective interaction. 

An alternative scenario was proposed, namely, the final 
state interaction of Y( 10860) decaying into Y(15\ 2S)7r + 7r~ 
can be realized via sub-processes Y( 10860) — > B^'B^ — > 
Y(1S ,2S)n + n~ l9T 4l"lll . It was also claimed that the anoma- 
lous Y(lS)n + 7T~ and Y(25')^ + 7r _ productions near the Y(5S) 
resonance receive reasonable explanations. The tetraquark 
interpretation presented in Ref. |]7|] is evidently not unique, 
and, moreover below we will show that the dipion invariant 
mass distribution and the angular distribution of y,(10890) — > 
Y(2S)n + 7r~ obtained with the teraquark picture proposed in 
Ref. 01 cannot explain the Belle data well. 

Using the formulas and parameters given in Ref. to fit 
the data given by the Belle Collaboration, we find obvious 
discrepancies. Namely, as shown in the following figure, em- 
ploying the formulas and parameters given in Ref. |7J|, we ob- 
tain the solid lines for the dipion invariant mass distribution 
and angular distribution for y,(10890) — > Y(2S )7r + 7r~ (see 
Fig. [TJ), and obviously they do not fit the data points which are 
marked in the figure. If, with the formulas given by the au- 
thors of Ref. 10], we fit the dipion invariant mass distribution 
(the dashed line) as the left-hand diagram of Fig. [TJto gain the 
model parameters, then applying the parameters to calculate 
the angular distribution dF/dcos 0, we would have the dashed 
line on the right-side of Fig. [TJ Inversely, if we first fit the 
differential width c/F/c/cos 9 to fix the parameters (the dotted- 
line on the right-hand side of Fig. [TJ, using those parameters 
would result in the dotted-line on the left-hand side diagram of 
Fig. [TJ All the results contradict the data; even the trend does 
not coincide. Therefore, it seems that the tetraquark scenario 
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does not work well to some extent. 
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FIG. 1 : (color online). Dipion invariant mass {m„ n ) distribution (left- 
hand side) and the cos 6 distribution (right-hand side) of Y(2S)7r + 7r~ 
production. The dots with error bars are the results measured by 
Belle. The solid lines denote the results reproduced with the param- 
eters shown in Table. II of Ref. [7J]. With the same formula as that 
in Ref. @|, we refit the experimental data. Here, the red dashed-line 
curves and blue dotted-line curves are the fitting results with param- 
eters {F = 0.933 ± 0.396, p = 0.692 ± 0.202, f a = 9.405 ± 2.409, 
4> = -0.460±0.245 Rad) and {F = 1.056± 1.348,/? = 0.467± 0.578, 
/„■ = 10.354± 15.163, (p = -1.785 ± 1.223 Rad), respectively, which 
correspond to the best fits to the dipion invariant mass distribution 
and the cos 9 distribution, respectively. 



In Ref. U, the authors considered the sub-processes 
T( 10860) fiWfiW _> r(lS,2S)n + n- which are sup- 
posed to be the final state interaction of T( 10860) — > 
Y(1S, 2S )n + 7i~ . They concluded that as the absorptive (imag- 
inary) part of the triangle diagrams dominate, one can expect 
an enhancement of 200 ~ 600 times compared to the partial 
widths of dipion emission Y(nS) — > Y(1S, 2S)n + n~ {n < 3). 
Moreover, even though for Y(4S), the BB channel is open, 
but due to the limit in phase space, the p-value suppresses the 
contribution from the BB intermediate states. The data show 
that the partial width of Y(4S) Y(1S, 2S)7r + 7T is only 2~ 4 
times larger than that of Y(3S) -> T(1S ,2S)n + n- by the Belle 
and Babar measurements. Thus it seems that the largeness of 
the dipion emission of Y( 10860) can be explained as coming 
from the final state interactions. However, in Ref. JH, the 
authors did not give a fit to the line shapes of the dipion in- 
variant mass distribution dr/dm n +„- and the differential width 
dF/d cos 9 explicitly, so they claimed that their results were 
roughly consistent with data. 

Even though the partial width of Y( 10860) 
T(\S,2S)n + n- is larger than that of T(nS) -> f{mS)n + n- 
(n = 3,4, in = 1,2) by two orders, there is still no com- 
pelling reason to ignore the contribution from the direct 
dipion transition process. Especially an interference between 
the direct transition and the contribution from the final 
state interaction (intermediate heavy-meson loops) may 
result in a sizable change to each scenario. We notice that 
r(Y(10860) -> Y(2S)7T + 7r) = 0.85 ± 0.07(stat) + 0.16(syst) 
MeV which is larger than T(Y( 10860) T{\S)n + ji-) = 
0.59 ± 0.04(stat) + 0.09(syst) MeV H, even though the later 
one has a larger final-state phase space. If the largeness is 
due to the composition of Y( 10860), this relation should be 
inverted. This intriguing experimental fact reported by Belle 
Jit] seems to be evidence to support an interference between 
the direct transition and the contribution from the final state 



interaction for Y( 10860) -> T{\S ,2S)n + n- . Obviously, the 
final state interaction that is realized via the heavy-meson 
loops is a simplified version of the complicated multichannel 
dynamics ifTHl . The direct dipion transition was dealt with in 
terms of the QCD multipole expansion method where there 
are two color-El transition and the two color fields eventually 
hadronize into two pions. Yan and Kuang |[l2tl established the 
theoretical framework for the multipole expansion method, 
where the intermediate state between the two El transitions 
the quark pair QQ resides in a color-octet where Q stands 
for heavy quark b or c. Here, we do not intend to calculate 
the contribution from the direct transition, but set it as an 
effective interaction with the free parameter, which will be 
fixed by fitting data. 

In this work, we suggest that the total decay amplitude of 
Y(55) — > Y(1S, 2S)n + n~ should include a few terms such as 
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where we take into account contributions from differ- 
ent intermediate resonances R to dipion, i.e., R = 
(cr(600), / (980), / 2 (1270)} to the Y(lS)7r + 7r channel and 
R = {cr(600), /o(980)} to the Y(2S)7r + 7T channel, which are 
allowed by the phase spaces. The phase angles <p^' are intro- 
duced, which will be fixed in this paper. 

In general, the decay amplitude of the direct production of 
T(5S) — > Y(15, 2S)n + n~ is expressed as 

M[Y(5S) ^ y(nS)( Pl )7T + (p 2 )7T-(p 3 )] 
= e^ nS) [\q 2 - 4AM) 2 (l + 7 ^f] 

VL \ q /JS-wave 

4(AM)^)(l-^)(cos^-I) 



which was first written by Novikov and Shifman in Ref. II 1 3fl 
while studying — > J/i[m + 7T~ decay, where the S-wave and 
D-wave contributions are distinguished by the subscripts S- 
wave and D-wave. AM denotes the mass difference between 
T(5S) and Y(nS). q 2 = (p 2 + pi) 1 = m^^- is the invariant 
mass of n + n~ . 6 is the angle between Y(55 ) and n~ in the n^ix" 
rest frame. In Ref. 10], Ali et al. also adopted the expression 
in Eq. (|2]) and introduced an extra form factor J?l = F/f 2 with 
f n = 130 MeV. 

As shown in Fig. [2 there are six diagrams correspond- 
ing to T(5S) decays into Y(nS)S and Y(«5')/2(1270) respec- 
tively, and S and /o(980) eventually turn into two pions; thus, 
as on-shell intermediate states they contribute to Y(5S) — > 
T{nS)n + ji~ . Such subsequent processes are attributed to the 
final state interaction. In this work, we adopt the effective La- 
grangian approach to write out the decay amplitudes for the 
diagrams in Fig. [2] where the relevant Lagrangians include 
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£rs-s = -igrra^d^dv&ifi +8djrf), (3) 
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+B*> 1 (r v d fi Bl t - 3 /i Y y S* vt )}, (4) 



and 



X^bwum = gsssSBB 1 ' - g s . s . 5 <SS*S t1 " (5) 



where B = (B°,B,B-) and {B 1 ) 7 = (B°,B + ,B+). Thus, 
the decay amplitudes corresponding to Figs. |2ja)-|2jf) are ex- 
pressed as 
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FIG. 2: The schematic diagrams for T(5S) decays into Y(«5)<S and 
T(«5)/ 2 (1270) (n = 1, 2) via B w meson loops. 



With above preparation, the amplitudes of Y(55') 
T{nS)n + ji~ via the re-scattering can be parameterized as 

3(*) d(*) 



M[T(5S) -> _» T(nS)(pik + (/7 2 K(^3)]s 

f (n) , (n) 1 ^(SS^HnS^s^Pi " P3 

y 1 q- - m s + inisis 

M[T(55) -» B ( *»B W -» T(«5)09 1 )7r + ( P2 )^(p 3 )] /2(1270) 
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corresponding to the contributions from the intermediate 
scalar states 5 = {cr(600),/o(980)} and the tensor meson 
/ 2 (1270). In the above equation, P^^iq) is defined as 

p%*</i) = \(r a ~g A " + 8 pB ~g Aa ) - \r A r B 

with g aB = g aB - q a qP lm\. Then the differential decay width 
reads as 
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TABLE I: The resonance parameters (in units of GeV) used in this 

work nun. 



Wt(5S) 


10.870 


m a 


0.526 rv 


0.302 


0.0010 


«T(IS) 


9.460 


m /o(980) 


0.980 r /0(980) 


0.070 




my(2S) 


10.024 


m h 


1.275 r f2 


0.185 


0.0005 



where mL = {p\ + po) 2 and m 2 nn = (p2 + ps) 2 - The factor 1 /3 
comes from an average over the polarizations of the initial 
Y(5S) state and in Ref. 0], this factor was missing. 

For the re-scattering process, the effective Lagrangian for 
coupling bottomonia to the bottomed mesons is determined 
based on the heavy quark effective theory ifTill . The coupling 
constants for Y(5S )B < -* ) B (t) are evaluated by fitting the partial 
decay widths while for T(nS)B {t) B ( *'> (n = 1, 2) and SB M B M , 
the coupling constants are directly taken from Ref. fl. In 
the re-scattering picture, for the Y(5S) — > Y(lS)7r + 7r~ pro- 
cess, the tensor meson /2(1270) should be included. Com- 
pared to the S -wave coupling f2B*B*, the D-wave couplings 
f2BB and f2BB* are negligible due to the high partial-wave 
suppression. The coupling between /i(1270) and B*B* has 
not been obtained from any measured reaction channels yet; 
thus, in present work, we treat this coupling constant as a free 
parameter to be fixed later. The coupling constants between 
the scalar mesons and the final n + n~ are g a „„ = 16.2 GeV~' 
and gf „„ - 2.40 GeV -1 , which are determined by fitting the 
corresponding partial widths. 

In our model, besides the phase angles tp { ^ between the re- 
scattering processes and the direct two-pion emission, just as 
indicated above, two more parameters J[ and k are introduced 
for accounting the contribution from the direct process as 
shown in Eq. ©. For calculating the re-scattering amplitudes, 
a form factor is employed to describe the off-shell effects of 
the exchanged mesons. In the calculations, the form factor 
takes the monopole form, i.e., t F(q 2 ) = (A 2 - m 2 E )l(q 2 - m 2 E ), 
where me is the mass of the exchanged meson in the 
fiWfiW -> f(nS)S, T(nS)f 2 transitions shown in Fig. |2] and 
A is usually reparameterized as A = niE + uAqcd- It is worth 
pointing out that such an adoption has a certain arbitrariness, 
but the value A, which manifests all the unknown information 
about the non-perturbative QCD effects and the inner struc- 
ture of the involved mesons, is determined by fitting data of 
various reactions; thus, it is believed that its value must fall in 
a reasonable range. Thus the arbitrariness is relatively allevi- 
ated. In present work, for T(5S) — > Y(nS)n + n~ (n — 1,2), 
a — 2 is adopted. The coefficients of the relevant Lorentz 
structures in Eqs. ([T2}-([T3) g™, g{*>, g™, (/ = 1, 2, 3) and 

g^Q are determined by calculating the hadronic loops. 

The m^n- and cos 9 distributions measured by the 
Belle Collaboration as well as the partial decay widths 
r T (5S)^T(is>r+tf- = 0.59 ± 0.04 ± 0.09 MeV and 
rT(55)^T(2S)^+^- = 0.89 + 0.07 ± 0.16 MeV, are taken as 
inputs to our work. All other input parameters, including the 
masses and widths of the involved particles, are listed in Table 
ID With the help of the MINUIT package, we fit the Belle 
data of Y(5S) — > Y(l S ,2S)tt + tt~ with the corresponding 
parameters being fixed and listed in Tables iHl and [TTT1 
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FIG. 3: (color online). Dipion invariant mass (m n * H -) distribution 
(left-hand side) and the cos 8 distribution (right-hand side) measured 
by Belle Q] f° r the final state T(1S)jt + jt~ (dots with error bars). The 
histograms are the best fit from our model. 



TABLE II: The parameters for Y(55) - 
by fitting the Belle data. Here = gf 2l 



Y(lS)/r + 7r that are gained 



Parameter 



Value 



F 

K 

8/2 



0.186 ±0.061 
0.459 ± 0.084 
12.361 ±20.109 



Parameter Value (Rad) 



^trt 600) 

,(1) 

<7q(980) 
'Tit 1270) 



-2.638 ± 0.735 
1.539 ±0.741 
-1.028 ±2.050 



The dipion invariant mass distribution dTldm n + K - and the 
angular distribution dT /d cos 9 measured by the Belle Collab- 
oration for Y(55') — > Y(lS)7r + 7r~ are shown in Fig. [3] The 
shaded histograms are the corresponding theoretical predic- 
tion by our model. The parameters for Y(55') — > Y(15')^ + 7r _ 
are listed in Table [ill yielding an integrated decay width of 
^r(5S)^r(\s)7T+7T- = 0.54 MeV. The consistency between our 
results and the Belle data indicates that our model can nat- 
urally describe the anomalous production rate of Y(15)7r + 7r" 
near the peak of Y(55' ) well, and, moreover, the predicted di- 
pion invariant mass distribution and the cos 9 distribution also 
coincide with the data. 

For Y(5S) — > Y(2S)7r + 7r~, we carry out a similar calcula- 
tion. The shaded histograms are our best fit to the experimen- 
tal data, and the corresponding parameters are listed in Table 
HIH The integrated decay width of r T(5S )^ T (2s)^+^- = 0.845 
MeV. Our results also confirm that the contribution from 
/o(980) is rather small compared to the Y(5S) -> Y(2S)7r + 7r 
process. 

However, one notices a discrepancy. As shown in Fig. [4] 
the dipion invariant mass distribution of the T(2S)n + 7r~ pro- 
duction near the peak of Y(55) is well reproduced by our 
model. However, applying the same fitting parameters, the 
predicted dF/d cos 9 of Y(55) — > Y(2S )7r + 7r~ (the histogram 
on the left panel of Fig. |4j, displays a different behavior from 
the Belle data for the Y(2S )ji + tt~ channel (dots with error bars 
in the right-hand diagram of Fig. |4). 

In summary, stimulated by the recent Belle observation of 



TABLE III: The fitted parameters for Y(55) -> Y(2S);r + 7r 



Parameter 


Value 


Parameter 


Value (Rad) 


F 

K 


2.315 ± 1.904 
0.572 ± 0.283 


tV(600) 
,(2) 

Y7 n (980) 


-0.297 ± 0.567 
-3.140 ±4.532 
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FIG. 4: (color online). The comparison between the fitting result 
(histogram) for T(55) -> Y(2S)7r + ;T and the Belle data (dots with 
error bars) 



anomalously large production rates of Y(1S, 2S)tt + ti~ near 
the peak of T( 10860), carefully studying the previous works 
along the line, we suggest that both the direct dipion emis- 
sion process and the processes via intermediate physical states 
which are the so-called final state interactions, contribute to 
the amplitude, and their interference results in the observed di- 
pion emission of Y(10860). In our scenario, the inverse rates 
r(Y(10860) -> t(2S)n + n-) > r(Y(10860) -> Y(lS)7r + 7r) 
can also be naturally understood, i.e., their rates are deter- 
mined by the interference between the contributions of the 
direct emission and the final interactions. 

Fitting the decay rates of Y(10860) -> r{\S, 2S)n + n~, 
we further theoretically investigate the dipion invariant mass 
and the cos<9 distributions of Y( 10860) -> Y(1S, 2S)n + n~ 
and make a comparison with the Belle data. Indeed, it is 
observed that if the final state interactions overwhelmingly 
dominate the transitions Y( 10860) -> Y(15) + tt+tT and 
Y(10860) -> Y(2S)7r + 7T, the lineshapes of the differential 
widths over the dipion invariant mass and cos 9 cannot be well 
fitted. It indicates that the interference effect plays a crucial 
role for fully understanding the Belle observation JU]. 

What is more important is that our model demonstrated 



in this paper shows that the tetraquark scenario proposed 
by Ali et al. I&d] does not provide a satisfactory under- 
standing of the anomalous e + e~ — » 4)(2S)tt + jt^ production at 
■\fs = 10.870 GeV. However, one cannot rule out that there 
might be a fraction of the tetraquark component in Y( 10860) 
that also contributes to the dipion transition. But, so far, it 
seems that a contribution from such an exotic state is not nec- 
essary for just understanding the Belle data. Instead, our study 
presented in this paper indicates that the Belle observation can 
be naturally explained by the interference between the direct 
dipion emission and the final state interactions. 

It is worth pointing out that in our model, our theo- 
retical predictions on both anomalous production rates of 
r(Y(10860) -> Y(2S);r + 7r) and T(Y(10860) -> Y(lS)7r + 7T) 
coincide with the data of the Belle collaboration, and also sat- 
isfactorily describe the dipion invariant mass and the cos 
distributions of Y( 10860) -> Y(15)tt + ^ as well as the dip- 
ion invariant mass distribution of Y( 10860) — > Y(2S)7r + ;r~. 
However, a new and intriguing puzzle is proposed since the 
predicted dT/dcos 6 of Y(10860) -> Y(2S)7r + 7T is inverse to 
the Belle data just presented in the right panels of Figs. [T]and 
[4] Associated with further theoretical exploration, future ex- 
perimental study from Belle-II and SuperB will be helpful to 
clarify this new puzzle and give a definite conclusion. 
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